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elliptic  partial  differential  equations. 


Numerical  predictions  of  the  backmlxed  flowfleld 
of  a methane-fired  laboratory  combustor  are  com- 
pared to  experimental  observations.  The  goal  Is 
to  clarify  the  mechanisms  responsible  for  pollutant 
formation  In  continuous  backmlxed  combustion  and  to 
Improve  predictive  methods.  The  numerical  des- 
cription adopts  a simplified  effective  viscosity 
model  and  a two-step  kinetic  reaction  scheme  for 
the  methane  oxidation.  A reaction  mechanism  to 
account  for  the  production  of  nitric  oxide  Is  also 
considered.  Laser  Doppler  veloclmetry  Is  used  to 
experimentally  measure  velocity  profiles  In  the 
flowfleld.  Conventional  gas  sampling  techniques 
are  used  to  measure  species  concentration  profiles. 

I.  Introduction 


Continuous  combustion  processes  represent  a 
category  of  combustion  characterized  by  an  un- 
interrupted, chemically  reacting  flow.  Continuous 
combustion  processes  feature  a complex  Interaction 
of  chemical  kinetics,  heat  transfer,  mass  diffusion 
and  fluid  mechanics.  An  Important  characteristic 
of  continuous  combustion  systems  Is  a zone  of  re- 
circulating flow  or  backmlxlng.  The  zone  serves 
as  a mechanism  for  flame  stabilization  by  back- 
mixing  hot  combustion  products  upstream  to  Ignite 
Incoming  fresh  reactants.  Gas  turbine  engines, 
utility  boilers  and  package  boilers  comprise  the 
largest  class  of  practical  continuous  combustion 
power  sources. 

Continuous  combustion  engineering  design 
practices  have  traditionally  been  based  on  consid- 
erations of  thermodynamic  efficiency  and  heat 
transfer  rates.  Present  operating  constraints, 
e.g.,  exhaust  composition  and  fuel  economy,  have 
necessitated  a more  detailed  study  of  the  processes 
occurring  In  the  combustion  system.  Because  the 
time-temperature  history  of  the  combustion  process 
plays  an  Important  role  In  determining  pollutant 
production  and  fuel  combustion  efficiency,  prior 
knowledge  of  concentration,  temperature  and 
velocity  distributions  Is  especially  valuable.  As 
a result,  the  design  and  operation  of  continuous 
coad>ustlon  systems  can  be  facilitated  by  predic- 
tions of  performance  by  way  of  an  analytical 
model. 

Theoretical  analyses  of  combustion-chamber  flow 
generally  incorporate  mathematical  models  of  the 
combustion  processes  In  a suitable  numerical  frame- 
work for  solving  the  governing  conservation 
equations.  Detailed  modeling  of  flows  dominated 


Progress  In  modeling  continuous  combustion 
systems  has  been  hampered  by  the  scarcity  of 
detailed  flowfleld  data  for  reacting,  recirculating 
flows.  Experimental  data  guide  the  development  of 
numerical  models  of  continuous  combustion  and  pro- 
vide a basis  for  their  evaluation.  A review  of 
the  state-of-the-art  In  continuous  combustor 
modeling  may  be  found  In  Osgerby  (1974).^ 

This  paper  presents  results  from  a continuing 
Investigation  designed  to  broaden  the  understanding 
of  the  physical  and  chemical  processes  occurring 
In  turbulent  combustion  representative  of  practical 
systems.  The  governing  equations  describing  a re- 
acting, turbulent  flowfleld  in  a backmlxed  methane/ 
air  combustor  are  solved  numerically.  Experimental 
data  are  obtained  and  compared  to  the  numerical 
description  of  the  flowfleld. 


II.  Background 

The  experimental  configuration  adopted  for  the 
present  study  is  a simplified  laboratory  burner 
consisting  of  a 2 In.  diameter  axlsysnetrlc  duct 
containing  an  aerodynamic  (opposed-jet)  flame- 
holder  as  shorn  in  Fig.  1.  The  opposed-jet  com- 
bustor (OJC)  provides  a versatile  bench-scale 
combustor  exhibiting  flame  characteristics  repre- 
sentative of  practical  combustion  systems.  The 
OJC  utilizes  a high  velocity  jet  stream  to  back- 
mix  the  primary  flow  of  reactants.**  The  flame 
is  stabllzed  near  the  stagnation  region  of  the 
main  and  jet  streams.  The  versatility  of  the 
opposed-jet  combustor  stems  from  the  controlling 
Influence  of  the  jet  composition  and  temperature 
on  the  reaction. 

The  numerical  procedure  used  to  solve  the  govern- 
ing equations  of  the  opposed-jet  flowfleld  Is 
based  upon  extended  versions  of  the  PISTEP  method 
of  Gosman  et  al  (1969)^.  The  numerical  procedure 
solves  simultaneously  the  basic  equations  describ- 
ing 2-D,  elliptic  flow  with  the  corresponding 
dependent  variables  listed  below: 


conservations  of  mass streamfunctlon, 

conservation  of  momentia vortlclty,^u>/r 

conservation  of  energy enthalpy,  h 

species  continuity mass  fraction,  m^ 

The  OJC  geometry  has  previously  been  Incorporated 
Into  the  computational  procedure  by  Samuelsen 
and  Starkman  (1972). ^ 


*Research  sponsored  In  part  by  the  Air  Force  Office 

of  Scientific  Research,  Air  Force  Systems  Command, 
USAF,  under  AFOSR  Grant  Number  2710-74 


**The  mass  flow  rate  of  the  jet  Is  much  less  than 
the  mass  flow  rate  of  the  main  stream 
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Turbulent  transport  properties  (exchange  co- 
efficients) are  evaluated  using  the  time-mean  or 
"effective"  quantities.  A simplified  effective 
viscosity  model  was  adopted  for  the  present 
investigation^ 

2/1  -V3  V3  ,0 

Meff  - KD  L p (mv^)  (1) 


recirculating  flowfield.  The  LDV  Is  operated  in  a 
forward  scatter  real^fringe  node  using  a SvW  He-Ne 
laser  light  source  and  a DISA  (Model  53L20)  fre- 
quency tracker  for  Doppler  signal  processing. 
Additional  features  of  the  LDV  system  have  been 
described  previously.^  Cold  flow  velocity  data 
are  used  as  a reference  base  for  evaluating  the 
system  aerodynamics,  the  numerical  predictive 
method,  and  the  associated  models  of  turbulence. 


Additional  exchange  coefficients  are  related  Co 

via  the  appropriate  Prandtl  or  Schmidt  number 

as  follows: 


The  effective  mass  exchange  model  ouiy  be  evalu- 
ated by  considering  a cold  flowfield  with  a non- 
uniform  concentration.  The  nonuniform  concen- 
tration is  established  by  introduction  of  a tracer 
into  the  jet.  SpAtlal  distributions  of  tracer 
concentrations  are  predicted  numerically  by  solving 
an  additional  species  conservation  equation  con- 
taining the  appropriate  mass  exchange  coefficient. 


Experimentally,  a carbon  monoxide  (CO)  tracer 
is  introduced  through  the  opposed-jet.  Local  CO 
concentrations  are  measured  by  NDIR  analysis  of 
gas  samples  extracted  by  conventional  probing. 
Comparison  of  the  predicted  and  experimentally 
determined  distributions  facilitate  identification 
of  deficiencies  In  the  coupled  mass  exchange  and 
effective  viscosity  models. 


Calculated  distributions  of  temperature  and 
species  mass  fraction  are  compared  to  experimental 
results  to  assess  the  coupled  models  of  turbulence 
and  kinetics  for  the  condition  of  reacting  flow. 

The  chemistry  model  incorporates  kinetic  mechanisms 
for  methane,  carbon  monoxide  and  nitric  oxide. 


A two-step  quasi-global  reaction  mechanism  for 
methane  oxidation^  is  adopted  for  the  initial  phase 
of  the  study: 

CH^  + I O2  — CO  + 2H2O  (3) 


Figure  1 Opposed-Jet  Conbustor 

a)  photograph 

b)  schematic  of  stream  function 
distribution 


Cold  Flow 


A detailed  description  of  the  numerical  formulation 
of  the  reaction  rate  expressions  may  be  found  in 
Samuelsen  and  Peck  (1972). ^ The  dependent  vari- 
ables for  the  hydrocarbon  system  Include  vortlclty. 
ii>/r,  streamf unction,  <)>,  methane  mass  fraction, 
m„„  , carbon  dioxide  mass  fraction,  m„^  , and 

COj 

enthalpy,  h.  Distributions  of  other  major  species 
H.O,  CO  and  0,  are  related  to  m.„  and  m via 

^ ^ UH.  LUa 

4 2 

elemental  mass  conservation. 


Consideration  of  cold  flow  conditions  enables 
a preliminary  evaluation  of  the  transport  models 
incorporated  Into  the  numerical  procedure.  A cold 
flow  solution  is  obtained  by  numerically  solving 
the  governing  equations  for  the  two  dependent 
variables  The  predicted  spatial  distributions 

of  the  time-average  (mean)  velocity  are  available 
for  comparison  to  the  experimental  measurements. 


Additional  calculations  were  conducted  to  pre- 
dict formation  of  the  pollutant  species,  nitric 
oxide  (NO) . Nitrogen  oxide  kinetics  were  based 
on  the  familiar  Zeldovlch^  mechanism: 

0 + Nj  — NO  + N (5) 

N + 0^  — ► NO  + 0 (6) 


Distributions  for  the  mean  velocity  are  measured  Noting  that  reaction  (5)  is  the  rate  limiting 

experimentally  by  using  a laser  Doppler  veloclmeter  step  simplifies  the  overall  reaction  rate  to: 

(LDV).  The  optical  dlagnoatic  technique  Is  desir- 
able because  It  does  not  disturb  the  characteristic 


Zkj  [Nj]  [o] 


Experimental  and  predicted  velocity  profiles  are 
presented  In  Fig,  2.  The  selected  profiles  empha- 
size the  flow  regime  dominated  by  recirculation. 


The  problem  remains  to  prescribe  the  O-atom  con- 
centration In  Eq.  (7).  Two  cases  are  considered 
In  the  present  analysis.  Equilibration  of  the 
elementary  reaction 


0 + 0 + M’‘  Oj  + M 


species  the  local  O-atom  concentration  via  the 
following  expression: 


Kg  [O^]** 


Fig.  2 Cold  Flow  Velocity  Profiles 
(•  exptl,-pred) 


Recent  evidence  indicates  that  the  assumption  of 
O/O2  equilibrium  underestimates  the  concentration 
of  monatomic  oxygen  atoms. The  so-called 
"super-equilibrium"  of  oxygen  atoms  was  incorporat- 
ed into  the  NO  kinetic  expression  following  the 
method  of  Iverach  et  al.  (1973).^  By  assuming 
that  the  reactions 


CO  + OH 

CO2  + H 

(10) 

and  H + 0 

- OH  +0 

(11) 

are  in  equilibrium,  the  following  expression  for 
the  concentration  of  oxygen  atoms  may  be  obtained: 


KipKii  [C0][0^] 


The  experimental  tests  conducted  to  complement 
the  nusierlcal  predictions  of  OJC  reacting  flow 
properties  utilized  stoichiometric  proportions  of 
premixed  methane/alr  as  reactants.  Temperature 
and  chemical  composition  measurements  were  obtain- 
ed at  the  OJC  exit  plane.  Temperature  measure- 
ments were  made  using  a Pt/Pt-13Z  Rh  thermocouple, 
Cas  samples  were  extracted  via  a water-cooled 
stainless  steel  probe  and  conveyed  through  a 
heated  teflon  sample  line  to  a packaged  exhaust 
gas  analysis  system.  The  Instrumentation  console 
provides  quantitative  analyses  of  CO  and  CO2 
(NDIR),  NO/NO  (chemiluminescence),  total  hydro- 
carbons (FID),  and  O2  (paramagnetic). 


IV.  Results  and  Discussion 


Cold  Flow 


The  cold  flow  studies  were  conducted  for  the 
following  test  conditions: 


Approach  velocity 
Jet  velocity 
Temperature 


• 25  fps 
■ 425  fps 


The  predicted  velocity  profiles  agree  qualita- 
tively with  experimental  results.  The  predicted 
location  of  the  stagnation  point  differs  sub- 
stantially from  the  location  Indicated  experimen- 
tally. The  Intense  mixing  processes  occurring  In 
the  recirculation  zone  suggest  employing  a refined 
model  of  turbulence  (eg.  two  parameter)  based  on 
the  turbulent  kinetic  energy  and  a prescribed 
mixing  length.  A thorough  evaluation  Is  currently 
limited  by  the  Inability  to  experimentally  obtain 
LDV  measurements  In  the  highly  turbulent  shear 
layer  between  the  jet  and  the  main  stream  or  In 
the  recirculation  zone  Itself. 

Results  of  the  tracer  studies  of  mass  transport 
In  a turbulent,  nonreacting  flow  are  shown  In 
Fig.  3.  Although  the  reported  trends  are  similar, 
the  calculated 'tracer  concentration  profiles  are 
quantitatively  higher  than  those  observed  experi- 
mentally. These  results  demonstrate  deficiencies 
In  the  mass  exchange  coefficient  and  more  Impor- 
tantly In  the  simplified  effective  viscosity  model 
used  In  the  current  study. 

a)  %C0  Predicted 


I 1 

r/x 


b)  %C0  ExperlmentRl 


Equivalence  Ratio,  " ♦j  * 


Fig.  3 Cold  Flow  Tracer  Concentration  Profiles 
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Results  of  the  reacting  flow  solution  for  the 
OJC  are  presented  In  Figs.  4 and  5.  The  conditions 
investigated  for  the  hot  flow  case  were  as  follows: 


Approach  velocity 
Jet  velocity 
Temperature  (Inlet) 
Equivalence  Ratio 


V ” 50  fps 
m 

Vj  = 425  fps 

1=1=  530°R 
m J 


b)  Taaperature.  *5: 


eet«:  radial  dlraetlon  expanded  2X 


Fig.  4 Hot  Flowfield  Property  Distributions 
(Predicted) 


Fig.  4 presents  profiles  of  stream  function, 
temperature  and  fuel  mass  fraction.  The  high 
temperature  reaction  zone  is  located  immediately 
downstream  of  the  stagnation  zone  and  parallels 
the  jet  wall  to  the  exit  plane  with  most  of  the 
fuel  being  consumed  prior  to  exhaust. 


a)  TfM  Caiien  Moaesld*  (00) 


Fig.  5 (continued) 


b)  PFN  Nitric  Oxide  (NO):  Caae  1 0/0,  tqlbra. 


c)  PPM  Nitric  Oxide  (NO):  Caae  2 O/O2  Noo-eqlbra. 


note:  radial  direction  expanded  2X 


Fig.  5 Hot  Flowfield  Property  Distributions 
(Predicted) 

Predicted  conc.intratlon  profiles  of  the  pollut- 
ant species  CO  and  NO  are  shown  In  Fig.  5.  Signif- 
icant CO  concentrations  persist  throughout  the 
flame  zone  and  reach  a maximum  In  the  transition 
(flame  front)  region  where  the  temperature  and 
oxygen  concentration  are  relatively  low.  The 
numerically  predicted  NO  profiles  for  Case  1 
indicate  that  the  maximum  concentrations  appear  in 
the  peak  temperature  zone  with  proportionate 
Increases  depending  on  residence  time.  The  non- 
equilibrium  condition  of  Case  2 results  in  a 30- 
fold  increase  in  NO  production. 


Fig,  6 OJC  Exhaust  HC  and  Temperature  Profiles 


Exhaust  plane  measurements  of  temperature  anu 
hydrocarbon  concentration  are  compared  to  predicted 
values  In  Fig.  6.  The  location  of  the  flame  front. 
Indicated  by  steep  temperature  and  concentration 
gradients,  coincides  with  theoretical  predictions. 
The  temperature  deviation  In  the  reaction  zone  may 
be  corrected  by  considering  radiation  losses  In 
both  the  numerical  and  experimental  results , 


Fig.  7 OJC  Exhaust  CO  and  NO  Profiles 

Fig.  7 presents  similar  results  for  CO  and 
Case  1 NO  concentrations.  The  Indicated  trends 
are  favorable,  but  the  magnitude  of  predicted  CO 
and  NO  levels  are  significantly  low.  Examining  the 
control] Ing  kinetic  rate  models  Indicates  that  the 
conversion  of  CO  ->  CO^  does  not  occur  as  rapidly 
as  Eq.  (4)  predicts.  The  observed  NO  levels  exceed 
those  predicted  via  the  Zeldovlch  mechanism  with 
0/0_  equilibrium,  yet  fall  well  below  the  non- 
equillbrium  calculations. 

In  general,  the  analytical  treatment  of  OJC 
reacting  flows  produced  results  simulating  gross 
flame  characteristics.  The  coupled  turbulent/ 
kinetic  models  adopted  for  the  initial  calculatlois 
require  further  refinements  before  adequate  experi- 
mental correlation  is  achieved. 


V.  Surmiary 

Numerical  predictions  of  the  turbulent,  back- 
mixed  flowfleld  of  a methane-fired  opposed-jet 
combustor  have  been  compared  to  experimental 
observations.  Although,  favorable  qualitative 
correlation  has  been  established,  deficiencies  In 
the  predictive  models  have  been  Identified.  A 
cold  flow  analysis  Identified  deficiencies  in  the 
simplified  turbulence  model  considered.  A hot 
flow  analysis  Identified  deficiencies  In  the 
coupled  turbulent/kinetic  models  adopted.  The 
results  from  this  preliminary  study  provide  a 
reference  base  for  further  experimental  verifica- 
tion and  formulation  of  refined  kinetic  and 
transport  mechanisms. 

VI.  References 

1.  Osgerby,  I.T.,  Literature  Review  of  Turbine 
Combustor  Modeling  and  Emissions,  AIAA  J,  1^ 
(6),  1974,  743-754. 

2.  Gosman,  A.D.,  Pun,  W.H.,  Runchal,  A.K., 
Spalding,  D.B.  and  Wolfshteln,  M,,  Heat  and 
Mass  Transfer  in  Recirculating  Flows,  Academic 
Press,  London,  1969. 


3.  Samuelsen,  G.S.  and  Starkman,  E.S.,  Analytical 
and  Experimental  Investigation  of  an  Ammonla/Alr 
Opposed  Reacting  Jet,  Comb.  Sci.  and  Techn. 

1972,  31-41. 

4.  Peck,  R.E.  and  Samuelsen,  G.S.,  Velocity  Measure- 
ments In  a Premixed  Flame  by  Laser  Anemometry, 
presented  at  the  Western  States  Section,  The 
Comb.  Inst.,  El  Segundo,  CA,  WSCI  73-28,  1973. 

5.  Williams,  G.C.,  Hottel,  H.C.,  and  Morgan,  A.C., 
The  Combustion  of  Methane  In  a Jet-Mixed  Reactor, 
12th  Symp.dntl.)  on  Comb.,  The  Comb.  Inst.,  1969. 

6.  Samuelsen,  G.S.  and  Peck,  R.E.,  Pollutant  Forma- 
tion In  Reacting  Flow  with  Recirculation, 
presented  at  the  Western  States  Section,  The 
Comb.  Inst.,  Monterey,  CA. , WSCI  72-21,  1972. 

7.  Zeldovlch,  Y.B.,  Sadovnlkov,  P.Y.  and  Frank 
Kamenetski,  D.A.,  The  Oxidation  of  Nitrogen  by 
Combustion,  Academy  of  Sciences  of  USSR, 

Institute  of  Chemical  Physics,  Moscow-Leningrad, 
1974. 

8.  Bowman,  C.T.,  Kinetics  of  Nitric  Oxide  Formation 
in  Combustion  Processes,  14th  Symp.  (Inti.)  on 
Comb.,  The  Combustion  Institute,  Pittsburg,  PA, 

1973,  729-738. 

9.  Iverach,  D.,  Basden,  K.S.  and  Kirov,  N.Y.,  For- 
mation of  Nitric  Oxide  in  Fuel-Lean  and  Fuel- 
Rich  Flames,  14th  Symp.  (Inti.)  on  Comb.,  The 
Combustion  Institute,  Pittsburg,  PA,  1973, 
767-775. 

10.  Saroflm,  A.F.  and  Pohl,  H.J.,  Kinetics  of 
Nitric  Oxide  Formation  in  Premixed  Laminar 
Flames,  14th  Symp.  (Inti.)  on  Comb.,  The 
Combustion  Institute,  Pittsburg,  PA,  1973, 
739-754. 

VII.  Nomenclature 

D combustor  diameter 

h stagnation  enthalpy 

K constant  in  effective  viscosity  expression 

Kj^  equilibrium  constant 

ki  kinetic  rate  constant 

L length  of  flowfleld 

m mass  flow  rate 

mj  mass  fraction  of  jth  species 

r radius 

T temperature 

V velocity 

r exchange  coefficient 

p viscosity 

p density 

Pj  Schmidt  Number 

Pj,  Prandtl  Number 

0 equivalence  ratio 

0 stream  function 

(V  vortlclty 

[i]  concentration  of  species  1 
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mechanisms  responsible  for  pollutant  formation  in  continuous  backmixed  combustio 
and  to  improve  predictive  methods.  The  numerical  description  adopts  a simplifie 
effective  viscosity  model  and  a two-step  kinetic  reaction  scheme  for  the  methane 
oxidation.  A reaction  mechanism  to  account  for  the  production  of  nitric  oxide 
is  also  considered.  Laser  Doppler  velocimetry  is  used  to  exerimental ly  measure 

velocity  profiles  in  the  flowfield.  Conventional  gas  sampling  techniques  are 
~ jto  measure  species  concentration  profilpR.,^  ^ 

DD  I JAN  73  1473 — EPiTiawar  Tfibv tnyo-BsoLETE — vf  vi 


'Pecles  concent • 
NOV  ii  li  obsolete 


UNCLASSIFIED 1 i . i 

TY  classification  of  this  page  flWicn  r>»i«  EnAfW)  V 

' 3 ^ ^ 


